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Abstract 

Wc present measurements of CP-violating parameters in B^{B^) fl'^s ^^'^ 
B^ rj'K^ decays based on a 41.8 fb~^ data sample collected at the T(4S') 
resonance with the Belle detector at the KEKB asymmetric-energy e~^e~ collider. 
We fully reconstruct one neutral B meson as a — ^ v'^s eigenstate and 
identify the flavor of the accompanying B from its decay products. From the distri- 
bution of proper time intervals between pairs of B meson decay points, we obtain 
the CP-violating asymmetry parameters = 0.28 ib 0.55(,stoi)lQ'Qg(,s7/.st), and 

A^,ko = 0.13 ± 0.32{stat)to^lisyst). We also reconstruct charged B^ rj'K^ 
decays and determine a direct-CP violating asymmetry value of Ar]'K± = (— 1-5 =b 
7.0{stat) ± 0.9{syst))%. 
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In the Kobayashi and Maskawa (KM) model, CP violation is incorporated as 
an irreducible complex phase in the weak-interaction quark mixing matrix [1]. 
Measurements of sin20i, where 0i = tt— arg(^^f^), from CP violation in 

cb 

S° — > ccK° decays by the Belle [2] and BaBar [3] collaborations established 
CP violation in the neutral B meson system that is consistent with KM ex- 
pectations. Measurements of sin20i based on other decay modes, especially 
charmless modes that are mediated by diagrams that contain virtual particle 
loops, provide important tests of the KM model. In this letter we describe 
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the first measurement of CP- violating asymmetries in the penguin-mediated 
decay v'^s^ ^iid an improved measurement of the direct CP- violating 

asymmetry in the decay B"^ — >■ tj'K^ [4]. 

The KM model predicts CP- violating asymmetries in the time-dependent 
rates for B^ and B^ decay to a common CP eigenstate, fcp- When an T(4S') 
decays into a B^B^ meson pair, the two mesons remain in a coherent p-wave 
state until one of them decays. The decay of one of the B mesons at time 
ttag to a final state, /tag, which distinguishes between P° and projects the 
accompanying B meson onto the opposite 6-flavor which decays to rj'Kg at 
time tcp- The decay rate has a time dependence given by 

pj^o(At) = -^:^{^ + Q ■ ^cp(At)}, (1) 

where Acp{At) is the time-dependent CP asymmetry, r^o is the B^ lifetime. 
At — tcp — ttag, and the 6-flavor q — 1) when the accompanying B 
meson is a B'^{B°). 

Within the framework of the Standard Model (SM), the charmless decay 
P° — > rj'Kg proceeds primarily via 6 — > s penguin diagrams; there is a small 
contribution from a color-suppressed b ^ u tree diagram, but that amplitude 
is expected to be only a few percent of that for the b ^ s penguin [5-7] . Thus, 
CP violation in this decay mode, to a good approximation, measures (f)i, and 
one can compare the result to the value measured for P° — > {cc)Kg. Phases 
from new physics in the penguin loop could show up as a difference between 
the two measured values [5,8]. Since the branching fraction of P — > rj'K ap- 
pears to be anomalously large [9-11], this mode is especially interesting to 
search for effects of an additional phase besides 4>i due to physics beyond the 
SM [5]. 

The time-dependent CP asymmetry can be expressed as 



r(po ^ 7]'Kl) - r(po ^ r]'K^s) 



r(po ^ r]'K°s) + r(po ^ r]'K^s) 
— '^v'K'^ cos(AmAt) -|- Sj^ij^o sin(AmAt), (2) 

where the CP-violating parameters Ajji^o and <S^'xg ^-^e given by 

'^^'^°S - |A|2 + 1' ^^'^°S - 1 + |A|2' ^'^^ 

in which A is a complex parameter that depends on both B^-B^ mixing and 
the decay amplitude for B^{B^) —>■ rj'Kg. The SM value for |A| is very nearly 
equal to the absolute value of the ratio of the P° to P° decay amphtudes. 
Therefore |A| 7^ 1, or equivalently Ari^x^ 0) indicates direct CP violation. 
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The direct CP asymmetry, A,fK±. in charged r]'K^ decays can also be 

investigated from the time-integrated decay rates of B~ versus B^. To a good 
approximation the CP- violating parameter, Sj^/xo, in the decay B^ — > rj'Kg 
is equal to the parameter sin20i, which can be directly compared with the 
value measured from S° — > {cc)K^ decays. 

The measurement reported here is based on a 41.8 fh^^ data sample, which 
contains 44.8 million BB pairs, collected with the Belle detector at the KEKB 
asymmetric energy (3.5 GeV on 8 GeV) e+e~ collider [12] operating at the 
T(45') resonance. At KEKB, the T(45') is produced with a Lorentz boost of 
= 0.425 nearly along the electron beam direction (z). Since the B^ and 
B^ mesons are approximately at rest in the T(4iS') center-of-mass system. At 
can be determined from the displacement in z between the fcp and /tag decay 
vertices: At ~ {zcp — Zta.gj/P'yc = Az/P'yc. 

The Belle detector is a large solid-angle general purpose magnetic spectrom- 
eter that consists of a three-concentric-layer silicon vertex detector (SVD), a 
50-layer central drift chamber (CDC), an array of aerogel threshold Cerenkov 
counters (ACC), time-of-flight scintillation counters (TOP), and an electro- 
magnetic calorimeter comprised of 8736 CsI(Tl) crystals (ECL) located inside 
a superconducting solenoid coil that provides a 1.5 T magnetic field. An iron 
flux-return located outside of the coil is instrumented to detect muons and 
Kl mesons (KLM). The detector is described in detail elsewhere [13]. 

The B — > r]'K event selection closely follows the method described in de- 
tail in a previously published report that describes the branching ratio mea- 
surement [11]; the time-dependent CP analysis is similar to that used for 
B^ — > TT~^7T~ and presented in Ref. [14]. 

For the analysis reported here, the event selection is shghtly modified for 
the time-dependent measurement from what described in Ref. [11] in order 
to retain more signal events and to keep better control of systematic errors. 
Candidate Kg tt+tt^ decays are reconstructed from pairs of oppositely 
charged tracks that are constrained to a common vertex and have an invariant 
mass that is within ±16 MeV/c^ of the nominal Kg mass. Two decay channels 
are used for r]' reconstruction: rj' — rjTi^TT~ (r)'^^^) with rj ^ •yy (rj-y^); and 
Tj' — >■ p^'-f {fj'p^) with p° TT^TT^. To increase the event yield, we omit the 
minimum transverse momentum pt requirement for charged tracks that was 
used in Ref. [11]. Instead, we require that all of the tracks have associated SVD 
hits and radial impact parameters \dr\ < 0.1 cm projected on the r-(j) plane. 
Better tracking requirement improves the vertex determination and gives less 
bias due to detector asymmetry. Particle identification information from the 
ACC, TOP and CDC dE/dx measurements are used to form a likelihood 
ratio in order to distinguish pions from kaons with at least 2.5(j separation 
for laboratory momenta up to 3.5 GeV/c. Candidate photons from r]-y-^ iVp-y) 
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decays are required to be isolated and have > 50 (100) MeV from the ECL 
measurement. The invariant mass of 77^^ candidates is required to be between 
500 MeV/c^ and 570 MeV/c^. A kinematic fit with an r] mass constraint is 
performed using the fitted vertex of the 7r"'"7r~ tracks from the rj' as the decay 
point. For 77^^ decays, the candidate p° mesons are reconstructed from pairs 
of vertex-constrained vr+vr^ tracks with an invariant mass between 550 and 
920 MeV/c^. The rj' candidates are required to have a reconstructed mass 
from 940 to 970 MeV/c^ for the tj'^^^ mode and 935 to 975 MeV/c^ for 7]'^ 
mode. Charged candidates are selected for the decay — > rj'K'^ based 
on the particle identification information described in Ref. [11]. 

Candidate B mesons are identified by combining rj' and Kg (K^) candidates to 
form the beam-constrained mass Mbc = \J (-^bcam)^ ~ (Pb"^)^) ^^'^ the energy 
difference AE = E'^^ — £'beam' where -Ebeam center of mass (cms) beam 

energy (nominally 5.29 GeV), and E^^ and are the cms energy and 
momentum of the B candidate. The Mbc signal region is defined as Mbc > 
5.27 GeV/c^; the AE signal region depends on the mode: it is —0.1 GeV 
< A£; < 0.08 GeV for r)'^^ and |A£;| < 0.06 GeV for ri'^^. 

We extract the signal yields with a simultaneous unbinned maximum- likelihood 
(ML) fit for both AE and Mbc. The signal distribution is a product of a Gaus- 
sian function in Mbc and a Gaussian plus a bifurcated Gaussian tail function in 
AE. The means and widths, as well as an overall normalization are the fitting 
parameters. The shapes of the background distributions, described below, are 
determined from sideband events in the 5.2 GeV/c^ < Mbc < 5.265 GeV/c^ 
and \AE\ < 0.25 GeV region, but with AE < -0.12 GeV or AE > 0.1 GeV 
for the rj'^^^ mode and \AE\ > 0.08 GeV for the rj'p^ mode; the difference is 
due to slight differences in resolution in the data as well as in the Monte Carlo 
simulation determined from the higher statistics r]'K^ mode. 

The dominant backgrounds are from e'^e~ — > qq continuum events (g = 
u,d,s,c). In this case, signal and background events can be partially sepa- 
rated by the event topology, which tends to be jet-like for qq continuum events 
and nearly isotropic for BB events. We use | cos 6*^1, the cosine of the angle 
between the thrust axis of the B candidate and that of the other particles 
in cms. The requirement | cos^rl < 0.9 rejects 50% of the background while 
retaining 90% of the signal events. This is sufficient for the 77^^^ mode, which 
is relatively clean. 

For the rj'p^ channel, an additional cut is applied using event shape variables. 
These variables include Sj_, which is the scalar sum of the transverse momenta 
of all particles outside a 45° cone around the candidate 1]' direction divided 
by the scalar sum of their total momenta, and five modified Fox- Wolfram 
moments [15], all combined into a single Fisher discriminant. In addition, we 
use cos 9b, the cosine of the angle between the B candidate fiight direction 
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Fig. 1. The AE and distributions for the rj'Kg candidates (histogram) and 
the ML fit results (sohd curve); (a) and (b) show the r/' — > rjinr mode, (c) and (d) 
the r]' — > /97 mode. The dashed curves are the results of the fit for qq continuum 
background; the dash-dot curves in (c) and (d) show the expected BB background 
in the mode from the Monte Carlo simulation. 



and the beam axis (z) in the cms, and a p helicity variable H, which is the 
cosine of the angle between the tt"*" momentum direction in the p rest frame 
and the p momentum direction in the rj' rest frame. All of these variables are 
combined to form a likelihood ratio LR = Ls/{Ls + Lqg), where Ls{qg) is the 
product of signal {qq) probability density functions. We determine Ls from 
Monte Carlo (MC) and Lqg from sideband data. We require LR > 0.5 for the 
7]'^ mode, which rejects 81% of the background and retains 80% of the signal. 

Sideband events are used to determine the shape of the continuum background 
distributions for the rj'^^^ mode and rj^^ modes separately. The AE shape is 
modelled by a linear background obtained from the Mbc sideband. The Mbc 
shape is modelled by the empirical function of Ref. [16]. From the ML fit we 
determine the total background fraction in the signal region to be 40.4% for 
the Tjl^^^Kg mode and 56.5% for the rj'p^Kg mode. The backgrounds are pre- 
dominantly qq continuum events with small contributions from other generic 
BB decays, mainly due to charm daughter particles from B ot B decay. These 
backgrounds, which are determined from a sample of MC simulated generic 
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B meson decays, are negligible for the rj'^^^K^ mode (smaller than 1%), but 
contribute 14% of the total background for the fj^^Kg mode. The larger back- 
ground for the Tj'p^Kg mode is due to wide width of p mass and combinatorial 
background of 7 and p to form rj' candidates. 

In the case of multiple candidates from the same event, we select the candidate 
with the best value from the 1]' mass constrained fit. The signal efficiency is 
determined from MC events to be 17.2 (16.5)% for rj'^^^K^ {r]'p^Kg). From the 
above- described ML fit, we obtain 27.7tl'l signal events in the rj'^^K'g mode 
and 45.5jl79 events in the v'py^s mode; the results of the fits are shown in 
Fig. 1. 

Leptons, A baryons, and charged pions and kaons that are not associated 
with the reconstructed 5° — > rj'K^ decay are used to identify the flavor of 
the accompanying B meson. We apply the same method used for the Belle 
sin 201 measurement [2]. We use two parameters, q and r, to represent the 
tagging information: g is a discrete variable that corresponds to the sign of 
the b quark charge and has the value q = +1 for a b (i.e., B^) tag, and q = —1 
for a b (5°) tag; r is an event-by-event, MC-determined flavor-tagging dilution 
factor that ranges from r — (no flavor discrimination) to r = 1 (unambiguous 
flavor assignment). The value of r is used to sort data into six intervals of r, 
according to flavor purity, the corresponding wrong tag fractions, wi{l = 1,6), 
and dilution factors (1 — 2wi) are determined for each r bin from data as 
described in Ref. [2]. 

The vertex positions for the rj'Kg and /tag decays are reconstructed using 
tracks that have at least one three-dimensional space point determined from 
associated r-(f) and z hits in the same SVD layer and one or more additional z 
hits in the other layers. Each vertex position is required to be consistent with 
the interaction point profile smeared in the r-(f) plane by the average trans- 
verse B meson decay length. The /tag vertex is determined from all remaining 
well reconstructed tracks after the B^ v'^s candidate tracks are removed. 
Tracks from other Kg candidates are not used. The MC simulation indicates 
that the typical combined track-finding and vertex-finding efficiency is 83% 
for r]'Kg decays and 95% for the /tag decays. The typical vertex resolution 
(rms) is 147 pm (89 pm) for the rj'^^^ {v'/ry) mode on the CP side, and 159 pui 
for the 2;tag on the tagging side. 

The proper-time interval resolution for the signal, i?sig(At), is obtained by 
convolving a sum of two Gaussians (a main component due to the SVD vertex 
resolution and charmed meson lifetimes, plus a tail component to account for 
poorly reconstructed tracks) with a function that takes into account the cms 
motion of the B mesons. The fraction in the main Gaussian is determined 
to be 0.97 ± 0.02 from a study of control samples of fully reconstructed B^ 
mesons [2]. The means (/^main, A*taii) and widths ((Tmain, crtaii) of the Gaussians 
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are calculated event-by-event from the rj'K^ and /tag vertex fit error matrices 
and the values from the fits [17]. The background resolution Rgg{At) has 
the same functional form but the parameters are obtained also event- by-event 
from the M^c and AE sideband events. We obtain the lifetimes for neutral 
and charged B mesons for the rj'K channels using the same procedure; the 
results [18] are consistent with the world average values. 

After vertexing, we find 77 candidate events with q = +1 (5°) fiavor tags and 
74 candidate events with g = -1 (W) for S° ri'K'^. Figures 2(a) and 2(b) 
show the observed At distribution for the two samples. 

We determine the CP- violating parameters, Ajjik° ^^'^ ^ri'Kl-i by performing 
an unbinned ML fit to the observed At distributions. For perfect resolution, 
the probability density function (pdf) for the signal, P^g, as a function of 
Aj^iK^ and iS^'xg is given by Eqs. 1 and 2 with q replaced by q{l — 2wi) to take 
into account the dilution due to mis-tagging. We fix r^o and Arrid at their 
world average values [19]. The pdf used for the qq background distribution is 

n,(At) = 2 |/r • -^2^ + (1 - fr) ■ S{At) I , (4) 

where is the background fraction with an effective hfetime Ug and 5 (At) is 
the Dirac delta function. We determine fr = 0.064 ± 0.022 and Ug = 2.24 ± 

0.37 ps from the sideband data. The pdf used to account for the small generic 
BB background for the rj'^ mode is V^-^i^At) = 6{At)/2. This background is 
mainly due to mis-reconstructed secondary (charm) decays, such as D^*^'^p'^, 
£)(*)Tp±^ /)(*)07r± and D^*Hv, etc. with D mesons decaying into a kaon and 
multiple pions or a p meson. 

We define the likelihood value for each event as: 

Pi = j {4nig(At', W;O^Sig(At. - At') + f[^gg{At')Rgg{Ati " At') 

+ /BB^BB(Ai')^BB(Aii - At')} dAt'. (5) 

Here fl {k = sig, qq or BB and / = 1, 6) are the weighted probability functions 
determined on an event-by-event basis as a function of AE and Mbc, properly 
normalized by the average signal and background fractions in the fitting region 
for each r interval for the signal rj'Kg, qq and BB events, defined as 



fk{AE, Mbc) = , pAE,Mbc , „l_W^E,Ml,o , „l _c^A_B,Mb, • (6) 
ysig-'^sig "T ilqq qq ^33 33 

Here gl^^ + g^^ + g^^-^ = 1 and i?^^'^bc q^.^ ^j^g shape functions for the signal 
{k — sig), continuum background {k = qq) and BB background {k — BB). 
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-8 -6 -4 -2 2 4 6 8 
At(ps) 

Fig. 2. The At distributions for the v'^s candidates in the signal region: (a) 

candidates with q = +1, i.e. the tag side is identified as B^; (b) candidates with 
q = —1; (c) r]'Kg yields after background subtraction; (d) the CP asymmetry for 
B^ Tj' after background subtraction; (e) the same asymmetry calculated for 
Tj' events. The curves in the figures show the results of the unbinned maximum 
likelihood fit. In (a) and (b), the solid, dashed and dash-dot curves are fit results for 
the total, qq background and BB background respectively. In (c) the solid (dashed) 
curve is for the q = +1 {q = —1) signal fit. In (d) and (e), the solid curve is the fit 
result for A^//^o and 5^//^o while the dashed curve corresponds to zero asymmetry. 



11 



The average event fractions, g^, are measured for rjirn and p7 separately in 
each r interval. 

In the fit, An'K^ and S^,j^o are free parameters that are determined by max- 
imizing the likelihood function C — Hi Pi, where the product is over all 
S° rj'K^ candidates. After maximizing the combined likelihood, the CP 
asymmetry parameters are determined from a total of 72.9 fj'Kg signal events 
(37.3 and 35.6 W-t&gs) to be: 

S^,KO = 0.28±0.55(stat)tl-^l(syst), 
A^,KO = 0.13±0.32{stat)t'ofe{syst). 

Figures 2(a) and (b) show the At distribution for B^- and i?°-tagged events 
together with the fit curves; the background-subtracted At distributions for 
signal-only are shown in Fig. 2(c). The errors on data points in Fig. 2(c) are 
statistical only and do not include the error associated with the subtracted 
background obtained by the fit. The background- subtracted time-dependent 
CP asymmetry between the B^- and i?°-tagged events is shown as a function 
of At in Fig. 2(d), with the result of the fit for S^ij^o and A^^k'^ superimposed. 

The systematic errors are summarized in Table 1 for S^,f^o and Aj^ip-o- The 
dominant sources for S^/^o are due to uncertainties in the signal and back- 
ground determination {AE and Mbc pdf's, event fractions and background 
shape), resolution functions, wrong tag fractions, vertexing, and the physics 
parameters (tbo and Am^). For Aj^^Kp the uncertainties in the signal and 
background determination and the vertexing are the leading components. We 
determine the systematic error due to physics parameters by repeating the fit 
after varying those parameters within their error taken from the world aver- 
age [19]. The systematic errors for wrong tag fractions, resolution functions, 
signal and background pdf's are estimated by repeating the fit after varying 
the parameters by ±1(7 determined from the data or MC. The systematic error 
due to vertexing is studied by varying the cut on vertex ot removing it and 
then repeat the fit. 

A number of checks are also performed. We analyze the rj'K"^ sample in the 
same way as rj'Kg. The raw asymmetry for the — > rj'K^ candidates is 
shown in Fig. 2(e). A fit to 230 candidate events yields S = 0.11 ± 0.29 and 
A = — 0.27±0.17, consistent with no asymmetry, as expected. We also examine 
the event yields and the At distributions for B^- and ^''-tagged events in the 
AE and Mbc sideband region. We find no significant difference between the 
two samples in both the rj'K'^ and rj'Kg modes. The average raw asymmetry 
from the sideband data is —0.009 ± 0.014, which is consistent with zero and 
indicates no bias. 

We also analyze the self-tagged rj'K^ event sample to search for a direct 
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AE (GeV) AE (GeV) 

Fig. 3. The AE projection plots for the rj'K^ candidates: (a) for rj' K~ , (b) for 
r]'K^. The cross hatched histograms indicate the ryvrvr channel; the solid histograms 
are the sum of both ij' decay channels (r/Trvr and p'j). The curves are the fitted 
backgrounds (dashed) and the sum of signal and background (solid). 

CP asymmetry. Here we do not use flavor-tag information from the other B 
meson. The rj'K^ sample is divided into ri'K~^ and ri'K~ samples. Since the 
CP asymmetry {Ar^'K± = {N{B-) - N{B+))/{N{B-) + N{B+)) for r]'K^ is 
time independent, we use the same event selection and analysis procedure as 
described in Ref. [11]. A simultaneous ML fit to Mbc and AE is performed 
for each sub-sample. The fitted numbers of signal events in the r/^^^ iVp-y) 

mode are 66.0t^°2° {7^-7-\ti) for 5" decays and 66.6^^^ i^T-^-wi) for B+ 
decays, respectively. The number of produced B~ and B~^ events are obtained 
by maximizing the product of the likelihoods for each submode, as shown in 
Fig. 3. From the fit we obtain ^^/x± = (—0.4 ± 11.0)% for the r]'^^^ mode and 
(— 2.5±10.0)% for rj'p^] here the errors are statistical only. Since the systematic 
errors on r]' reconstruction and the number of BB events cancel in the ratio, 
the systematic uncertainty in Ari'K^ comes mainly from the reconstruction 
efficiency of charged kaons and the ML fit. The asymmetry in the efficiency 
is studied using inclusive charged kaons in the same kinematic range as the 
signal. The uncertainty due to fitting is measured by varying the parameters 
of the fit functions. We find the systematic errors in Arj'K^ ^'I'e 0.85% from 
reconstruction and 0.47% from the ML fit. The combined Ar^'K± for the 
B^ — > fj'K'^ decay is determined to be 

A^,K± = (-1-5 ± 7.0{stat) ± 0.9{syst))%, 

which is consistent with zero. Combining the errors in quadrature and as- 
suming a Gaussian distribution, we find the 90% confidence level interval is 
—0.13 < Ari'K^ < 0.10, which is a factor of two more restrictive than our 
previous measurement [11]. 

In summary, we have measured the CP asymmetry parameters in B^[B^) 
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rj'Kg and rj'K^ decays based on a 41.8 fb~^ data sample collected with 

the Belle detector. The result for decay-time-integrated direct CP asymme- 
try, Ar^'K± = (—1-5 ± 7.0{stat) ± 0.9{syst))%, is small and consistent with 
zero. Our results for the time-dependent asymmetry parameters S^/j^o = 

0.28 ± 0.55{stat)tomisyst) and A^^ko = 0.13±0.32{stat)to°ol{syst) are the 
first measurements of CP asymmetry parameters related to (pi with a charm- 
less decay mode. In the SM, to a good approximation the value of 
should be equal to sin 20i measured in B^ — > ccK^ decays, where the current 
world average is 0.78 ± 0.08 [20]. With a much larger data sample, we will 
significantly reduce the uncertainty in S^ji^o and impose tight constraints on 
phases from new physics beyond the Standard Model. 
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Table 1 



Summary of systematic 


errors. 




Source 






T^o and Aijid 


+0.029 
-0.034 


+0.006 
-0.003 


Wrong tag fractions 


-0.028 
-0.025 


+0.013 
-0.012 


Resolution functions 


+0.037 
-0.033 


+0.011 
-0.008 


Vertexing 


+0.018 
-0.033 


+0.062 
-0.031 


Event fractions 


+0.029 
-0.024 


+0.012 
-0.010 


Background shapes 


+0.013 
-0.017 


+0.004 
-0.003 


AE and Mbc pdf's 


+0.032 
-0.028 


+0.063 
-0.054 


Total 


+0.07 
-0.08 


+0.09 
-0.06 
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